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tion  of  a model  which  describes  the  generation  and  dispersal  of  aerosols 


resulting  from  explosive  and  other  smoke  munitions.  Here  we  have  extracted 
the  essential  formalism  and  organized  It  Into  a step-by-step  procedure 
which  will  permit  the  computation  of  (1)  the  spatial  history  of  the  aerosol 
cloud,  (2)  aerosol  concentration,  and  (3)  smoke  (obscuring  aerosol)  trans- 


mittance. In  view  of  the  exothermic  properties  of  white/red  phosphorus 
(WP/RP)  the  clouds  produced  by  WP/RP  are  treated  separately;  that  Is  the 
effect  of  buoyancy  Is  Included. 

We  have  listed  here  some  pertinent  data.  Whenever  extensive  data 
are  required  the  reader  Is  referred  to  Volume  I,  e.g.,  list  of  munition 
characteristics. 

With  respect  to  the  optical  properties  of  smoke  materials  these 
are  still  In  a state  of  flux  and,  therefore  should  be  treated  accordingly. 

The  adaption  of  this  model  to  combat  simulation  models  is  not  treated 
here;  however,  the  necessary  Input  data  are  listed.  Further,  we  hav  not 
altered  the  original  (July  1977)  suggestion  that  the  smoke  cloud  may  be 
treated  as  a transient  loss  of  Intervisibility,  see  the  discussion  in 
Volume  I. 

In  Section  2 we  consider  the  methods  to  compute  the  concentration 
of  a smoke  cloud  and  In  Section  3 the  transmission  and  other  optical 
effects  for  both  passive  and  active  systems  are  considered. 

In  Appendix  A,  we  give  a discussion  of  the  asymptotic  values  for 
plume  heights. 
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2.  COMPUTATION  SCHEME  FOR  CLOUD  CONCENTRATION 

In  this  section  a procedure  Is  outlined  to  permit  the  determina- 
tion of  the  concentration,  C (gm/m*),  path  length  In  the  cloud,  L (m), 
and,  the  important  parameter,  C*L  (gm/m2)*.  The  units  used  here  are 
defined;  however,  we  utilize  primarily  grams,  meters  and  seconds.  We 
Initiate  the  discussion  with  the  requisite  Input  data  and  then  proceed 
to  the  procedure  proper. 

2.1  Input  Data 

A.  Scenario 

Local  target  and  observer  locations 
Munition  deployment  plan 
Munition  rate 

Sensor  type  and  operational  wavelength  band 

B.  Smoke  Munitions 

Type  (Mk  No,  Caliber,  etc.) 

Fill  weight  (gm) 

Fill  material  (WP,  RP,  HC,  Oil,  etc.) 

Burn  time  and  rate 

C.  Meteorological  Conditions 
Time  of  day  ^ 

Cloud  cover  J To  establish  stability  category 

Mean  wind  speed  (m/sec)  - surface  to  10  m. 

Temperature,  T(°C)  - 0 10  m and  0.5  m above  ground  level 
Wind  direction 
Relative  humidity  (RH  %) 

D.  Optical  Properties 

Scattering,  absorption  and  extinction  coefficients 
Visibility,  If  eq.  (28)  and  Tables  9 and  10  are  used. 
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2.2  Determine  Concentration  for  Single  or  Separate  Multi -Munitions 

2.2.1  Establish  Stability  Category 

(a)  Find  AT  » T(10  m)  - T( .5  m) 

If  AT  =•  0.098°C  then  CAT  D. 

AT>0.098°C  enter  Table  1 with  Input  data  and  find 
stability  category 

Table  1 


PASQUILL  STABILITY  CATEGORIES 


Surface  wind 
speed  at  10m 
(m/sec) 


Stronc 

A 

A-B 


Day 

Insolation 
Moderate  SI i ght 


Night  thinly 
overcast 
or  > 4/8 
low  cloud 


i 3/8  cloud 


NOTE:  If  potential  temperature,  0(T)  * AT  + r,  is  zero  we  have  a neutral 
(D)  category*  1 .e. , 

0 (T)  =•  0 = -0.098  + 0.098 


2.2.2  Establish  Cloud  Dimensions 

(a)  Non-Exothermic  (HC,  oil,  etc.) 

x * u-t,  u given  In  Input  data 

^ | with  stability  category  enter  Table  2 and  find 
2 * appropriate  expressions.  For  more  precise  value 
of  plume  heights  see  Appendix  G,  Vol  T! 

where  x * coordinate  parallel  to  u (m) 
u ■ time  mean  wlndspeed  (m/sec) 
t ■ time  (sec) 
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l 

i 

Table  2 

>* 

1 

I 

PLUME  DIMENSIONS 

♦ 

Category  Width,  y (m)  Height,  z (m) 

Remarks  * 

i T 

A 

9.1  + .419X  2.73  + .137X 

Sunny  day 

«% 

B 

9.1  + .328X  2.73  + .llx 

Day,  broken  clouds 

it 

C 

9.1  + .237X  2.73  + .073X 

Overcast  day/night 

u 

D 

9.1  + .2X  2.73+.06X 

Neutral 

n 

u 

E 

9.1  + .18X  2.73+.055X 

Evening/early  morning 

F 

9.1  + .146X  2.73+.046X 

Evening/early  morning 

I i 

u 

I r ? 

- 

For  CAT  D (see  Appendix  a) 

7 1 50  n > n 

zmax  * — ’ 0 

1 

I 

0 

For  CAT  E,  F (see  Appendix  a) 

'max  - =>° 

0 

(b)  Exothermic  (WP/RP) 

n 

Compute  x,  y,  z as  above 

1 

To  calculate  the  effect  of  heat  on  plume  rl 

se.  Ah,  l.e.,  to  find 

[i 

Z * z + Ah,  we  follow  the  following  procedures: 

l u 

0 

D 

Categories  A,  B & C: 

Ah  * 1.6  F^ 

Q x 

f 

F * 3.7  x 10‘5H 

0 

H ■ 800x  gms  of  WP/RP  fill  wt  (gms)  f burn  time  In  secs* 

Z|nax  occurs  when  C = Q(gms  of  smoke)/V(cloud  volume  in  m3) 

n 

a 0.1  g/m3 

1 

Li 

NOTE:  will  occur  for  times  less  than  40  secs.  Method 

1 

D 

to  calculate  V Is  given  later 

D 

fi 

* Burn  time  for  explosive  devices  Is  1 sec.  Other  munitions  must  be 
obtained  from  PM- Smoke  or  CSL. 

L 

Neutral  Category  D: 

Ah  a LlEl  x*2/-  [.4  + .64(f*)  + 2.2(|*)2][1  + .8  fj 
x*  ■ 0.52F*^hs*®;  (hs  * 1 meter) 

(Ah)max  occurs  for  3x*  s x s 5x* 


Stable  Categories  E,  F & G: 

use  Ah  for  unstable  category  except  limit  to  a maximum  value 
(AlOmax  * 2.9  (F/G5)1/3 

xfflax  * 2.4  u/S*5  (see  Appendix  A) 

where  5 » I 

T az 

g * 9.8  m/sec2 

T » mean  ambient  temperature  (°K) 

0 * potential  temperature 


[ 

0 

L! 

0 


0 


2.2.3  Establish  Mass  of  Smoke  Produced 

(a)  Specify  smoke  type  and  relative  humidity 
enter  appropriate  equation  below  and  find  yield 

Y (WP ) * .003  (RH  - 40)1,67  + 2.9 

Y(ZnCl2)  ■ .051  (RH  - 50) *85  + 1 

Y (A1C1 3/H2SO4)  - .016  (RH  - 20)1,24  + 1.4 

Y(Fog  011)  » 1 

(b)  Mass  of  smoke  * Yield  x Munition  Fill  Wt  (gm) 
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2.2.4  Smoke  Concentration 


(a)  Concentration  for  HC,  Fog  Oil,  Metal  Chlorides,  etc 


Mass  of  smoke 


denominator  ■ volume  of  semi -cone 


(b)  Concentration  for  WP/RP 


Mass  of  smoke 
i x(t)  y(t)  Z(t) 


denominator  3 volume  of  elllpsoid/4. 

NOTE:  Establish  value  of  C(t)  * 0.1  gm/m3,  thereafter  Z(t)  remains 

constant.  Reiterate  on  C(t).  Note:  x and  y continue  to  Increase 


2.3  Smoke  Path  Length,  L 

(a)  For  a single  munition;  draw  line-of-sight  and  find  length 
through  smoke  cloud. 


(b)  For  a multi -munition  (mortars,  artillery  shells,  bombs,  smoke 
Same  as  (a). 


2.4  Calculate  C*L 


(a)  For  single  munition,  use  appropriate  result  from  (2.2)  and 


(b)  For  multi -munition,  use  appropriate  result  (for  single  event) 

n 

from  (2.2)  and  (2.3)  and  then  sum  the  result  from  n-rounds  Z CX»LX 


For  oblique  viewing,  the  llne-of-sight  may  initially  Intercept 
only  a portion  of  the  expanding  smoke  cloud. 


3. 


CONCENTRATION  FOR  CONTINUOUS  LINEAR  SMOKE  SOURCE  OF  LENGTH,  l 


I 


(a)  Stability  category  - Find  y(t);  z or  Z. 

(b)  Calculate  concentration 

„ Mass  of  Smoke  1 

**t#U  ^ .zU)-(Hk) 


(g/m3) 


Mass  of  Smoke  ■ Total  weight  x yield  (gms);  2.2.3. 


I - length  of  smoke  (m) 
t ■ burn  time  of  smoke  material  (sec) 
Q = mean  wind  speed  (m/sec) 

Z(t)  = cloud  height  (m) 

k ■ cloud  width  increment;  y(t)/£ 


(c)  Smoke  Path  Length,  L 

Draw  line-of-sight  through  smoke  cloud. 

(d)  Concentration  x Path  Length 
(2)  x (3)  = C-L 


3.1  Two  dimensional  or  Several  Linear  Smoke  Sources 
See  Volume  I,  pps  73-76  and  Appendix  D. 


0 

D 

0 

0 

0 


4.  TRANSMISSION  AND  OPTICAL  EFFECTS 

We  consider  here  only  the  transmission  through  smoke  and  other 
optical  effects,  e.g.,  glare,  produced  by  a smoke  cloud.  The  remainder 
of  the  path  between  target  and  observer  is  not  treated  here  since 
computer  programs  exist  to  perform  this  calculation.  We  will  first 
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consider  passive  systems  and,  then  attempt  to  outline  a procedure  for 
active  (laser)  systems. 


I 


I 


I 


l 


0 


n 


4.1 


Transmission,  T:  Passive  Systems 

(a)  Determine  sensor  wavelength  or  passband-lookup  extinction, 
scattering  and  absorption  coefficients  for  smoke  type. 

(b)  Establish  sensor  field-of-view  and  target  radiation  angle. 


NOTES:  (1)  For  viewing  systems,  e.g.,  FLIR,  use  instantaneous 
field-of-view  and  not  total  field-of-view. 


(2)  Passive  targets  are  normally  LAMBERTIAN.  Searchlights 
etc.  are  collimated  sources. 


4.1.1  LAMBERTIAN  Radiator 

Transmission  = (l+as*L*D)  exp-[ae(m)*C*L] 

The  value  for  D is  given  below  for  several  wavelengths  and  three 


The  correction  for  scattering  is  only  important  for 
visible  and  near  infrared  systems. 

Whenever  the  coefficients  vary  within  a passband,  then 
integration  over  the  band  should  be  performed. 


vi 

receiver  fields-of-view,  <t>,  = 

• .5,  5 & JjO. 

X(um) 

TIP5  5*  TCF' 

. . 

0.5 

.40  .60 

1.0 

.30  .50 

5.0 

.10  .20 

0 

10.0 

.05  .10 

n 

NOTE:  D = correction 

for  scattering 

IJ 

ct$  * scattering 

coefficient  (m_1) 

0 

<*e(m)  * extinction 

o 

coefficient  (m  /gm) 

4.1.2 


Collimated  Sources 

A simpler  form  for  transmittance  may  be  used  for  collimated 


I 

I 

I 

! i 


sources 


T 3 [1  + Os»L*<J>]  • exp[-ae*L] 


where  <}>  3 receiver  f.o.v.  in  degrees/ 180. 


NOTE:  A laser  transmitter  Is  a collimated  source  and  the 
spectral  extinction  coefficient  ae(X)  * ae(X,m)*C 
in  units  of  m"1 


NOTICE:  Establish  system  criterion  for  transmission  as  to  lower  permissible 

limit,  i.e.,  if  T < limit  then  terminate  calculation. 


5. 

CONTRAST,  C 

n 

5.1 

Visual  and  Near  Infrared:  Daytime 

We  use  the  Blackwell  definition  of  contrast,  namely 

. 

0 

iTarget  Brightness  - Background  Brightness! 
background  Brightness 

0 

C is  the  contrast  at  the  scene.  The  contrast  Cr,  at  observer 

0 

Is  defined  by  the  expression 

cr*co1'4  (e®6*-!)]’1 

0 
r ] 

R » distance  of  cloud  edge  to  observer  (m) 

Apply  test  if  Cr < 5% 

l j 

(a)  For  sunny  day  S«/G  = .2/8 

(b)  Other  conditions  S^/G  3 1/6 

0 

n 

Value  of  8 is  given  In  Table  3 for  different  values 
of  as/cte,  cig'L. 

1 

0 

|0 

9 

L 

IMP 

I 


Table  3 
VALUES  OF  0 


I 

VL 

T 

ois/oie  0 

0.5 

1.0  2 

3 

l 

.1  0 

.00 

.0  .015 

.02 

.3 

.02 

.03  .045 

.05 

•ft 

.5 

.04 

.07  .09 

.10 

.7 

.08 

.12  .17 

.19 

M 

.9 

.14 

.25  .33 

.37 

gy 

.95 

.22 

.35  .43 

.47 

- 

1.0 

.22 

.35  .50 

.59 

Li 

The  values  for  as/ae  for  HC  and 

WP  smoke  clouds 

are  tabulated  below: 

0 

Table  4 

WAVELENGTH 

DEPENDENCE  OF  cts/oie 

[1 

(J 

X(ym) 

V“e 

HC 

WP 

n 

< 1 

1 

1 

Li 

2.0 

1 

- 

2.5 

.974 

- 

0 

3.0 

.474 

.249 

3.5 

.841 

.242 

o 

4.0 

.894 

.274 

4.5 

.876 

.255 

ri 

5.0 

.835 

.256 

Vi 

8.0 

.485 

.042 

n 

8.5 

.392 

.043 

LJ 

9.0 

.309 

.050 

n 

9.5 

.245 

.036 

li 

10.0 

.184 

.067 

B 

10 

A. 


Table  4 (Cont'd) 


as/ae 


A(pm) 

hC 

w r- 

10.5 

.129 

.126 

11.0 

.088 

.123 

11.5 

.055 

.133 

12.0 

.033 

.102 

12.5 

.029 

- 

13.0 

.023 

- 

NOTE:  If  ae*L>3,  transmission  <_  5%  and  therefore  no  need  to  calculate 
contrast  reduction. 

If  Cr  > 5%  then  utilize  simple  result  or,  the  one  described  in 
Volume  1,  pps  53-55. 

5.2  Near  Infrared:  Nighttime 

Use  simple  procedure  with  night  sky  spectral  irradiance 
curve  in  place  of  solar  irradiance. 

5.3  Mid-IR  (3- 5pm) 

We  describe  here  only  the  glare  contribution  due  to  solar 
scattering  from  the  smoke  cloud  as  the  contrast  loss  due  to  the  ambient 
environment  is  treated  elsewhere.  We  recommend  here  the  use  of  the 
expression 


H(S)  ■ solar  Irradiance 

P(0)  » differential  scattering/unit  cloud  area/steradlan 

cosi.e  ■ cosine  of  angle  of  Incidence  and  exltance 

(x)  ■ atmospheric  spectral  transmittance  from 
source  to  cloud 

Nx(e)  ■ spectral  steradiance  exiting  from  cloud 
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NOTES:  1.  For  passband,  integration  of  all  parameters  on  right  hand 
side  is  necessary. 

2.  Values  of  P(©),  etc,  must  be  obtained  from  Mie  program  or 
experiment. 

5.4  LWIR  (8-1 3ym) ; (3-5ym,  Nighttime) 

Cloud  glare  is  not  important,  no  calculation  Is  required. 

6.  LASER  SYSTEMS 

6.1  Transmission 

See  4.1  for  one  way  pass;  i.e.,  cloud  between  transmitter  and 

receiver. 


When  cloud  exists  between  receiver  and  laser  illuminated  target 
use  Section  4.1.1.  Add  solar  background  (cloud  backscatter)  clutter 
where  the  receiver  pass-band  (AX)  must  be  used. 

6.2  Backscatter 

See  Section  5.3  for  appropriate  expression.  The  value  for 
P(©)  must  be  measured  or  computed  from  Mie  scattering  theory.  The  latter 
is  available  from  C.S.L.  Aberdeen  (J.  Vervier). 


7.  THERMAL  EFFECTS 

In  the  above  described  model  we  have  assumed  the  cloud  concen- 
tration to  be  homogeneous.  Estimates  of  the  material  and  vertical  dis- 
tributions are  discussed  In  Volume  I;  however,  these  distributions  may 
not  present  sufficiently  severe  gradients  to  affect  thermal  viewing  or 
laser  systems.  A detailed  analysis  of  the  cloud  structure  will  be 
performed  in  an  upgraded  model. 
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A problem  to  thermal  or  laser  systems  that  can  be  readily 
identified  is  the  multiple  source;  that  Is,  if  a number  of  phosphorus 
sources  exist  false  targets  will  be  generated  in  the  field  of  view. 
This  effect  Is  Identical  to  scene  clutter  and  should  be  treated 
accordingly. 


Appendix  A 


Plume  Rise  Limiting  Values 


The  purpose  of  this  appendix  is  to  clarify  the  several  rules 
governing  the  maximum  values  that  the  exothermic  smoke  plume  can  achieve 
and  explain  the  additional  rules  added  here  in  Section  2.2.2  for  non- 
exothermic  smokes/aerosols.  We  will  first  consider  the  non-exothermic 
obscurant. 

In  Table  2,  Section  2.2.2,  we  define  the  plume  rise  as  a monotonic 
increasing  linear  function  of  x,  i.e.,  2 = A + Bx.  For  atmospheric 
conditions  defined  as  Pasquill  Category  D (neutral)  and  Categories  E & 

F (stable)  we  have  established,  on  the  basis  of  data  analysis,  that 
the  cloud  stabilizes  at  specific  heights;  that  is,  the  maximum  height 
can  be  defined  at  each  value  of  x,  the  horizontal  displacement,  by  A/u 
where  A is  a constant.  A set  of  values  is  given  in  Table  A.l  below. 

TABLE  A.l 

MAXIMUM  VALUES  FOR  PLUME  RISE 

u.zmax  (m2/s) 

Stability  Category 


x(m) 

Neutral 

E&F 

250 

60 

25 

500 

105 

55 

1000 

150 

75 

A-l 

Let  us  now  consider  the  case  of  a neutral  category  with  u 3 5 m/s. 


I Then,  from  Table  A.l  the  zmax  at  x 3 250rn  Is  12m  whereas  from  Table  2, 

z * 2.73  + .06x  * 2.73  ♦ 15  * 17.7.  If  we  attribute  the  value  of  2.73 
to  source  dispersion  at  x-0  then  the  difference  Is  -25%.  One  finds  by 

O 

the  same  procedure  that  the  difference  at  100  m Is  '15%.  On  the  other 
hand  If  we  let  x * 1000  m we  have  from  Table  A.l  that  z,^  3 30  m and 
from  Table  2,z  ■ 60  m or,  100%  difference.  We  assert  that  the  plume 
rise  equations  In  Table  2 represent  the  linear  portion  of  the  rise  and 
as  the  plume  becomes  bent  due  to  the  horizontal  wind  (mean  value  u) 
other  expressions,  as  given  In  Table  A.l,  are  applicable.  Hence  It  Is 
recommended  that  Table  A.l  be  used  for  x^250  m. 

For  exothermic  smokes,  we  stated  that  the  maximum  rise  due  to 
buoyancy  under  stable  atmospheric  conditions  Is 

v 

Ah  3 2.9  (F/us)  , (1) 


where  F Is  the  rate  of  heat  released,  and  should  occur  at  a maximum 
- -k 

value  of  x 3 2.4  u/s  . The  general  rise  has  been  defined  by 

h 3 1.6  V/3  xZ/3(F1/s)  (2) 

1/ 

where  we  added  the  extra  multiplier  F , i.e.,  we  altered  the  exponent 
of  F from  Va  to  l/2,  to  agree  with  experimental  results  for  WP  munitions. 
If  we  take  the  conventional  value  (F^3)  in  eq.  (2)  and  equate  eq.  (1), 
we  find  that  the  maximum  rise  will  occr  at 

Xmax  ‘ 2.4u/sl/2  (3) 

which  Is  independent  of  "F".  On  the  other  hand,  the  explicit  use  of 
eq.  (2)  gives  for  x|||ax 

Li  xiMX  • 2.45/s1'4  • FV"  (4) 

0 

r 


A- 2 


w 

1 1 1 

II 

The  appearance  of  "F"  In  eq.  (3)  raises  some  questions.  First,  as 
j F-+0,  (and  Ah-KJ)  a condition  which  probably  is  not  physically 

meaningful.  Secondly,  when  F-*-»,  x^^  we  have  a condition  that  is 
questionable;  that  Is,  the  buoyancy  will  dissipate  at  x~0.  The  reader 
Is  reminded  that  the  established  criterion,  namely  the  plume  rise  due 
to  buoyancy  be  terminated  when  the  smoke  concentration  drops  to  0.1  g/m3; 
•*  is  also  applicable. 

In  the  next  edition  of  this  model,  the  questions  raised  here 
relative  to  complete  internal  self-consistency  will  be  addressed  and 

resolved. 
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